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Abstract. - We study the steady-state current in a minimal model for a quantum dot dominated 
by charge fluctuations and analytically describe the time evolution into this state. The current is 
driven by a finite bias voltage V across the dot, and two different renormalization group methods 
are used to treat small to intermediate local Coulomb interactions. The corresponding flow equa- 
tions can be solved analytically which allows to identify all microscopic cutoff scales. Exploring 
the entire parameter space we find rich non-equilibrium physics which cannot be understood by 
simply considering the bias voltage as an infrared cutoff. For the experimentally relevant case of 
left-right asymmetric couplings, the current generically shows a power-law suppression for large 
V. The relaxation dynamics towards the steady state features characteristic oscillations as well 
as an interplay of exponential and power-law decay. 
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Introduction. — Recent progress in the ability to en- 
gineer nanostructured devices has opened new possibili- 
ties for studying the finite-bias transport characteristics 
of such systems. As the electrons occupying the nanos- 
tructure are spatially confined, local Coulomb correla- 
tions strongly affect the physics, and understanding non- 
equilibrium phenomena in systems with local two-particle 
interactions is therefore of fundamental importance. In 
an attempt to investigate simplified cases first, one can 
distinguish between situations in which either charge or 
spin fluctuations dominate. The latter case is described 
by the Kondo model, and progress in understanding its 
non-equilibrium physics was made recently (for a review 
see Ref. [T]). We here consider the other situation and 
study a minimal model for a quantum dot dominated by 
charge fluctuations — the interacting resonant level model 
(IRLM). It describes a spinless localized level at energy 
e coupled to two leads by electron hoppings t a and local 
Coulomb repulsions u a (see Fig.[T]). The lead electrons are 
assumed to be (effectively) non-interacting and held at two 
different chemical potentials /i Q = ±V/2, with a = L, R 
denoting the left and right lead and V being the bias volt- 
age. 



The steady-state current I of the IRLM was studied 
intensively during the last few years using various tech- 
niques including the scattering Bcthe Ansatz [5] , perturba- 
tive and numerical renormalization group (NRG) methods 
[3J , the Hcrshfield F-operator [3] , the time-dependent den- 
sity matrix renormalization group (tDMRG) method [5] as 
well as sophisticated field theory approaches 016]. These 
studies were mostly performed at the non-generic point 
of particle-hole and left-right symmetry, which can hardly 
be realized in experiments. It was concluded that at suffi- 
ciently large V, I decreases as a power law. Similar power 
laws were found in equilibrium and it was suggested that V 
is just another infrared energy cutoff (in addition to, e.g., 
t a or temperature), leading to the speculation that the 
IRLM does not contain any interesting non-equilibrium 
physics [3J . Exploring the entire parameter space we show 
analytically that this conclusion is too restrictive. We un- 
cover rich non-equilibrium physics beyond the situation 
where the voltage V acts as a simple low-energy cutoff as- 
sociated with a power-law behavior of the current. How- 
ever, in the limit of strong left-right asymmetry, which can 
be easily realized experimentally, we find generic power- 
law scaling of / for large V, in particular also away from 
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Fig. 1: The interacting resonant level model discussed in this 
work. 



particle-hole symmetry. 

In addition, we provide an analytic description of the 
relaxation dynamics of the system into the steady state 
after switching on the level-lead coupling at time t = 
0. Two different relaxation rates control the exponen- 
tial decay which is accompanied by oscillatory behavior 
with a voltage-dependent frequency and power-law de- 
cay with an exponent depending on u a . Describing the 
time evolution of a locally correlated electron system is as 
challenging as understanding the non-equilibrium steady 
state current. Various numerical techniques like time- 
dependent NRG and tDMRG g], an iterative path- 
integral method [3], and a non-equilibrium Monte Carlo 
approach [TU] were developed. Certain exactly solvable 
models were discussed [11J . and a perturbative renormal- 
ization group (RG) method [THUS] as well as a flow equa- 
tion approach [T3] were applied. However, these studies 
do not cover charge-fluctuating, correlated quantum dots. 

In this Letter, we use two RG methods to investigate 
the IRLM. While both are bound to the case of weak 
Coulomb interactions, they are complementary in other 
aspects. Within the functional RG (FRG), which was re- 
cently extended to non-equilibrium [T5], the steady state 
can be studied for arbitrary system parameters. In par- 
ticular, this allows for a comparison to highly accurate 
tDMRG data obtained for hoppings which are too large 
to be deep in the scaling limit [5], the latter being real- 
ized for large band width and small t a . For small inter- 
actions, we find excellent agreement (see Fig. EJa)). In 
the scaling limit, the FRG results and the ones obtained 
by the real-time renormalization group in frequency space 
(RTRG-FS) [J coincide (see Figs. [2] and !}. The latter 
method was earlier applied to systems dominated by spin 
fluctuations [H]. In contrast to the FRG, RTRG-FS can 
only be used in the scaling limit, but on the other hand 
allows for an analytical description not only of the steady 
state but also of the relaxation dynamics. The combined 
use of both RG approaches leads to a reliable and com- 
prehensive picture of the non-equilibrium physics under 
consideration. In particular, we identify the various mi- 
croscopic cutoff scales, which is essential for the precise 
determination of the scaling behavior of observables. 

Model and RG equations. — The Hamiltonian of 
the IRLM (sec Fig. [TJ) is given by H = Hi + H d + H c , 
where Hi = J2ka( e k + ^a)a\ a aka describes two semi- 
infinite fermionic leads which are held at Pl/r. = 



respectively. Standard second quantized notation is used, 
and the energies e& a are restricted to a finite band of width 
B. In the scaling limit, the details of the frequency depen- 
dence of the lead local density of states p a do not play any 
role as long as it is sufficiently regular for energies of the 
order of V and smaller. When comparing to tDMRG data 
[5], we employ the semi-circular p a {oj) associated with 
simple tight-binding chains (which are used in tDMRG). 
The dot Hamiltonian reads Hd = en with n — c^c, 
and this single fermionic level is coupled to the leads via 



H r . 
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where : . . . : denotes normal-ordering. We stress that in 
contrast to other studies, the coupling to the leads is al- 
lowed to be asymmetric, which is the situation generically 
expected in experiments. Furthermore, we do not only 
focus on the particle- hole symmetric point e = 0. 

Within both RG approaches, coupled differential equa- 
tions for the flow of the effective system parameters as a 
function of an infrared cutoff A holding up to leading order 
in u a can be derived. Aiming at an analytic discussion, 
it is instructive to consider simplified flow equations for 
the renormalized steady-state rates r a whose bare values 
are given by F Q a = 2ixp a t 2 a . In the scaling limit, both RG 
approaches give the same functional form 



dA 



~2U a r D 



A + r/2 



(p a - e)2 + (A + r/2)= 



(1) 



with U a = p a u a being the dimensionless interaction, and 
r = ^ Q r Q . The renormalization of the level position e 
is small and will be neglected. The RG flow (JXJ) is cut off 
at the scale A c = max{|^ Q — e|, T/2}, and an approximate 
solution for r Q is given by 



A„ 
A, 



2U a 



(2) 



where Ao ~ B denotes the initial cutoff. At large voltages 
V 3> r, wc distinguish between the off-resonance \V — 
2e| > r and the on-resonance V — 2e situation (peak in 
conductance; see Fig. H]). In the latter case, the relevant 
energy scales cutting off the flow are T/2 for Fl and V for 
Tr. 

Similiar to the Kondo model the cutoff parameter is the 
maximum of the distance to the resonance and the cor- 
responding decay rate, i.e. in our case max(|e ± V/2|,r). 
There is, however, an important difference. For the Kondo 
model, even at resonance V = h (the latter being the 
magnetic field), there is a weak-coupling expansion pa- 
rameter, namely the dimensionless exchange coupling cut 
off at max(V, h). For the IRLM, at e = ±V/2, the tunnel- 
ing is not a weak-coupling expansion parameter since F is 
not dimensionless. This fact constitutes an essential dif- 
ference between the description of resonance phenomena 
in models with charge and spin fluctuations. 

The full FRG and RTRG-FS flow equations are pre- 
sented in Refs. [T7| and [IB], respectively and can easily 
be solved on a computer. If not mentioned otherwise, 
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Fig. 2: (Color online) I(V) for the symmetric model with Ul = 
Ur = 0.1/-7T obtained from the numerical solution of the full 
RTRG-FS (lines) and FRG (symbols) equations; crosses: e = 0; 
stars: e = V/2. 

the results shown in the Figures were obtained in this way 
(for a comparison to the analytic solution of the simplified 
equations see Fig. 03b)). 

In the scaling limit where Ao — > oo and T° — > 0, the 
dependence on bare parameters vanishes, and all quan- 
tities can be expressed in terms of the invariant scale 
t k = E a T%, with T% = T a (2A /T K ) 2U « , and the asym- 
metry parameter c 2 = T K /Tg. Thus, at V = 



T = T h 
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which has the solution T = Tk in the symmetric case 
(Ul = Ur, c = 1). The corresponding equation for T at 
finite V in the off- (on-) resonance situation is obtained 
by replacing T's on the right-hand side of ^ by V — 2e 
and V + 2e (by T and 2V). As a result, the rates r a are 
gencrically characterized by power laws with interaction- 
dependent exponents. 



Steady-state quantities. 

the stationary state reads 



{U) = 9 + Z 
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and the static susceptibility is defined as x = — ^r|e=o- 
In the symmetric case and at V = 0, one obtains Xv=o = 
2/(7rr) [T5], which can be used to define the physical scale 
Tk = ^(Xy=o) _1 even away from the scaling limit. The 
stationary current can directly be computed from the rates 
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For V 3> r and off resonance, this expression simplifies to 
I ~ r L r H /r and thus 
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Fig. 3: (Color online) a) Comparison of FRG results (lines) 
and tDMRG data [5] (symbols) of I(V) for u L = u R = 0.3B/4 
and t L =t R = 0.5B/4 at e = 0. b) RTRG-FS results for I(V) 
for Ul/r = (l±7)0.1/7r, and e = 0; analytic result |(2]| inserted 
in (J5]) (dashed lines) mostly hidden by solution of the full flow 
equations (solid lines); 7 = 0.75, 0.5, 0.25, (corresponding to 
c 2 = 21.4, 7.9, 2.8, 1) from top to bottom. 



Whereas for Ul — Ur = U and V » e (e.g., at the 
particle-hole symmetric point e = 0) the current is al- 
ways governed by a power law I(V) oc V~ 2U in agreement 
with earlier studies [U[5], this does not hold in general for 
asymmetric Coulomb interactions generically realized in 
experiments. In this case the two terms in the denomi- 
nator of © arc typically of the same order of magnitude. 
Only if in addition to Ul 7^ Ur the asymmetry in the 
bare rates is large (c C 1 or c > 1), the power-law be- 
havior of /(V) is recovered (with exponents 2Ul or 2Ur, 
respectively). In the on-resonance case e — V/2 where the 
conductance G = dl/dV has a maximum (see Fig. 2]), the 
current is given by 

(Zs.) 2C/i (Ik) 2C/B 



TlTr = Tk 
2T 2 



'IT 



1 (Tj^^Ur 1 
c \ 2V ) 



+ c 2 



(7) 



c ^ \V+2n\ J 



In contrast to the off-resonance situation, / does not follow 
a power law even in the left-right symmetric model (see 

( 4 ) Fig. E]). Only for very large V (or for c 3> 1), the second 
— o- term in the denominator of ([7]) can be neglected and / oc 

V~ 2Ur [2D]. Thus, the voltage V cannot be interpreted as 
a simple infrared cutoff both for e = ±V/2 and Ul ^ Ur 
and the physics in non-equilibrium is far more complex 
than in the linear-response limit |21| . 

The analytic results ©, (JB]), and Q derived from ap- 
proximate FRG and RTRG-FS flow equations arc con- 

(5) firmed by solving the full RG equations numerically. The 
current for the left-right symmetric model exhibits a 
power-law decay / oc V~ 2U and thus a constant logarith- 
mic derivative only in the off- resonance case (see Fig. [2]). 
Fig. EJb) illustrates for e = and different coupling asym- 
metries that the current from the full RTRG-FS flow equa- 
ls) tion is captured by the analytic solution for the rates ([2]) 

inserted in ([5]). Moreover, the FRG compares nicely with 
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Fig. 4: (Color online) Conductance G(e) = dl/dV (lines: 
RTRG-FS, symbols: FRG) in the symmetric model with Ul = 
Ur = 0.1/vr. 



accurate tDMRG reference results obtained for large hop- 
pings (see Fig.[3ja)). Another transport property of exper- 
imental interest is the conductance G, which as a function 
of the gate voltage e most importantly features the men- 
tioned resonance at e — ±V/2 as the voltage becomes large 
(see Fig.H|). As before, both RG frameworks give agreeing 
numerical results for arbitrary values of V/Tk and e/Tx, 
thus altogether providing reliable tools to study quantum 
dot systems out of equilibrium. 

Time evolution. - The RTRG-FS allows for study- 
ing the time evolution towards the steady state. To this 
end, we initially prepare the system in a state described 
by p(t < 0) = ptjj PlPr, where pffl is an arbitrary initial 
density matrix of the dot and Pl/r are grandcanonical 
distributions of the leads. At time t = 0, the coupling H c 
is suddenly switched on and transient dynamics of pjj sets 
in. The latter can be fully described in terms of T a as a 
function of a Laplace variable z which has to be incorpo- 
rated |18j . By analytically solving an approximation to 
these RG equations, one can obtain closed integral repre- 
sentations both for the dot occupation and the cur- 
rent I(t) by virtue of inverse Laplace transform |13j . Nu- 
merical results for the time evolution are shown in Fig. [5] 
We restrict ourselves to the left-right symmetric model for 
simplicity. The long-time behavior away from resonance 
(i.e., at e, V, \e - V/2\ > T K , l/t) is given by 
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where (h) follows from and Ti « T as well as 

T 2 ~ Ti/2 are two decay rates. Whereas Ti describes 
the charge relaxation process, T 2 is the broadening of the 
local level e induced by the coupling to the leads, i.e. it 
describes the relaxation of nondiagonal elements of the lo- 
cal density matrix with respect to the charge states. We 
note that the dephasing rate = T 2 — T\/2 ~ 0(U) is 
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Fig. 5: (Color online) Time evolution of the dot occupation 
(n(t)) and the current I(t) for U L = Ur = 0.1/tt, e = 10T K , 
and the initial condition (n(0)) = [22]. At times t < 2/T K 
we observe oscillating behavior. 



due to pure potential fluctuations on the dot and increases 
for large Coulomb interactions. Most notable character- 
istics of the time evolution of both (n(t)) as well as the 
current I(t) are that (i) the relaxation towards the station- 
ary value is governed by both decay rates, (ii) the voltage 
appears as an important energy scale for the dynamics set- 
ting the frequency of an oscillatory behavior, and (hi) the 
exponential decay is accompanied by an algebraic decay 
cx t 2U ^ 1 . The last result is of particular importance for 
applications in error correction schemes of quantum infor- 
mation processing as it contrasts the standard assumption 
of a purely exponential decay [23] ■ We also note that in 
the short-time dynamics a reversal of the current can oc- 
cur (see the dashed-dotted curve in Fig. [5]Jb)). This effect 
is due to very strong charge fluctuations in the transient 
state, thus being impossible in systems with spin or orbital 
fluctuations |13j . Another interesting observation is that 
in the resonance case (dashed line) current oscillations are 
fully damped. 

Conclusion. — We have studied non-equilibrium 
transport properties of a spinless single-level quantum dot 
coupled to leads via tunneling and Coulomb interaction, 
representing a fundamental model to describe the effects 
of charge fluctuations. Using two different RG methods 
we have presented analytic results in the entire param- 
eter regime and concluded that the steady-state current 
I(V) exhibits a power law only in specific cases. The 
one of highest experimental relevance is the situation of 
strong asymmetries in the tunneling couplings, where we 
generically observed a power law for large bias voltages V. 
Furthermore, the time evolution towards the steady state 
was studied. We found exponential decay on two different 
scales accompanied by voltage-dependent oscillations and 
power laws with interaction-dependent exponents. 
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